clear

set mem 500m

spgrid, shape(sq) res(w1) xrange(0 100) yrange(0 100) cells("payoutcells.dta") points("payoutgrid.dta") replace

use "payoutgrid.dta", clear

drop spgrid_id spgrid_status spgrid_xcoord spgrid_ycoord
rename spgrid_xdim x
rename spgrid_ydim y


*******************************************
*** Generate all necessary variables    ***
*******************************************

	* IDEALPOINTS
		generate s1_opt_x1=22
		generate s1_opt_y1=35
		generate s1_opt_x2=17
		generate s1_opt_y2=45
		generate s1_opt_x3=49
		generate s1_opt_y3=51
		generate s1_opt_x4=82
		generate s1_opt_y4=45
		generate s1_opt_x5=80
		generate s1_opt_y5=36


	* Distance from every point in grid to ideal point

		foreach b in ///
				1 2 3 4 5 ///
				{
				generate s1_dist`b' = ((x-s1_opt_x`b')^2+(y-s1_opt_y`b')^2)^0.5
				}
			

	* Computation of resulting payouts (in points or tokens)

		// Treatment 1
		// Player 1, Treatment 1
		generate t1_points1 = 1000 - 5 * s1_dist1
		replace t1_points1 = 9320*exp(-1* (s1_dist1 / 14.8)) if s1_dist1 > 36

		// Player 2, Treatment 1
		generate t1_points2 = 1000 - 5 * s1_dist2
		replace t1_points2 = 9320*exp(-1* (s1_dist2 / 14.8)) if s1_dist2 > 36

		// Player 3, Treatment 1
		generate t1_points3 = 1000 - 5 * s1_dist3
		replace t1_points3 = 9320*exp(-1* (s1_dist3 / 14.8)) if s1_dist3 > 36

		// Player 4, Treatment 1
		generate t1_points4 = 1000 - 5 * s1_dist4
		replace t1_points4 = 9320*exp(-1* (s1_dist4 / 14.8)) if s1_dist4 > 36

		// Player 5, Treatment 1
		generate t1_points5 = 1000 - 5 * s1_dist5
		replace t1_points5 = 9320*exp(-1* (s1_dist5 / 14.8)) if s1_dist5 > 36


		// Treatment 2
		// Player 1, Treatment 2
		generate t2_points1 = 1000 - 5 * s1_dist1
		replace t2_points1 = 22000*exp(-1* (s1_dist1 / 14.8)) if s1_dist1 > 50

		// Player 2, Treatment 2
		generate t2_points2 = 1000 - 5 * s1_dist2
		replace t2_points2 = 29800*exp(-1* (s1_dist2 / 14.8)) if s1_dist2 > 55

		// Player 3, Treatment 2
		generate t2_points3 = 1000 - 5 * s1_dist3
		replace t2_points3 = 9320*exp(-1* (s1_dist3 / 14.8)) if s1_dist3 > 36

		// Player 4, Treatment 2
		generate t2_points4 = 1000 - 5 * s1_dist4
		replace t2_points4 = 2110*exp(-1* (s1_dist4 / 14.8)) if s1_dist4 > 12

		// Player 5, Treatment 2
		generate t2_points5 = 1000 - 5 * s1_dist5
		replace t2_points5 = 1865*exp(-1* (s1_dist5 / 14.8)) if s1_dist5 > 10


		// Treatment 3
		// Player 1, Treatment 3
		generate t3_points1 = 1000 - 5 * s1_dist1
		replace t3_points1 = 1865*exp(-1* (s1_dist1 / 14.8)) if s1_dist1 > 10

		// Player 2, Treatment 3
		generate t3_points2 = 1000 - 5 * s1_dist2
		replace t3_points2 = 2110*exp(-1* (s1_dist2 / 14.8)) if s1_dist2 > 12

		// Player 3, Treatment 3
		generate t3_points3 = 1000 - 5 * s1_dist3
		replace t3_points3 = 9320*exp(-1* (s1_dist3 / 14.8)) if s1_dist3 > 36

		// Player 4, Treatment 3
		generate t3_points4 = 1000 - 5 * s1_dist4
		replace t3_points4 = 29800*exp(-1* (s1_dist4 / 14.8)) if s1_dist4 > 55

		// Player 5, Treatment 3
		generate t3_points5 = 1000 - 5 * s1_dist5
		replace t3_points5 = 22000*exp(-1* (s1_dist5 / 14.8)) if s1_dist5 > 50



		// Total points
		foreach a in ///
			1 2 3 ///
			{
				generate t`a'tot_points = (t`a'_points1 + t`a'_points2 + t`a'_points3 + t`a'_points4 + t`a'_points5)
				}
	
		// Social optimum "Efficient point"
		foreach a in ///
			1 2 3 ///
			{
				egen max_t`a'tot_points1 = max(t`a'tot_points)
				generate max_t`a'tot_points = 0
				replace max_t`a'tot_points = 1 if max_t`a'tot_points1==t`a'tot_points
				drop max_t`a'tot_points1
				}

		// Average earnings
		foreach a in ///
			1 2 3 ///
			{
				generate t`a'avg_points = (t`a'tot_points)/5
				}

		// Inequality: Standard deviation of payouts
		foreach a in ///
			1 2 3 ///
			{
				gen t`a'_sd_points = ///
					(((t`a'_points1 - t`a'avg_points)^2 ///
					+ (t`a'_points2 - t`a'avg_points)^2 ///
					+ (t`a'_points3 - t`a'avg_points)^2 ///
					+ (t`a'_points4 - t`a'avg_points)^2 ///
					+ (t`a'_points5 - t`a'avg_points)^2)/5)^0.5
		}

		// Maximin
		foreach a in ///
			1 2 3 ///
			{
				egen t`a'_min = rowmin(t`a'_points1 t`a'_points2 t`a'_points3 t`a'_points4 t`a'_points5)
				egen max_t`a'_min1 = max(t`a'_min)
				generate max_t`a'_min = 0
				replace max_t`a'_min = 1 if max_t`a'_min1==t`a'_min
				drop max_t`a'_min1
				}
	
		// Maximum of individual payouts
		foreach a in ///
			1 2 3 ///
			{
				egen t`a'_max = rowmax(t`a'_points1 t`a'_points2 t`a'_points3 t`a'_points4 t`a'_points5)
				}


		foreach x in A B C D E A* B* C* D* E*{
		capture gen player_`x'="`x'"
		}
		
		// Coordinates inside UC-SET
		gen ucset = 0
		replace ucset = 1 ///
			if ///
			(y == 51 & x == 49) | ///
			(y == 50 & x >= 44 & x<=54) | ///
			(y == 49 & x >= 39 & x<=59.5) | ///
			(y == 48 & x >= 37.5 & x<=60.3) | ///
			(y == 47 & x >= 37.5 & x<=60.8) | ///
			(y == 46 & x >= 37.5 & x<=60.8) | ///
			(y == 45 & x >= 37.5 & x<=60.8) | ///
			(y == 44 & x >= 37.5 & x<=60.2) | ///
			(y == 43 & x >= 38 & x<=60.8) | ///
			(y == 42 & x >= 38 & x<=60.8) | ///
			(y == 41 & x >= 38 & x<=60.8) | ///
			(y == 40 & x >= 38 & x<=59.8) | ///
			(y == 39 & x >= 39 & x<=59.3) | ///
			(y == 38 & x >= 39 & x<=59.3) | ///
			(y == 37 & x >= 40 & x<=58.3) | ///
			(y == 36 & x >= 41 & x<=56.8) | ///
			(y == 35 & x >= 42 & x<=55)


			
********************************************************************************
*** PLOTS
********************************************************************************	

	*** Figure 2: Configuration of ideal points and location of the uncovered set
		
		twoway ///
			(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
			(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
			(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
			(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
			(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
			(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
			(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
			(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
			(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
			(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
			(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
			(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
			(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
			(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
			(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
			(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
			(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
			(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
			(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
			(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
			(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
			(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
			(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
			(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
			(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
			(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
			(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
			(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
			(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
			(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
			(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
			(pcarrowi 35 22 45 17, lcolor(gs8) lpattern(dash) msize(zero)) ///
			(pcarrowi 45 17 51 49, lcolor(gs8) lpattern(dash) msize(zero)) ///
			(pcarrowi 51 49 45 82, lcolor(gs8) lpattern(dash) msize(zero)) ///
			(pcarrowi 45 82 36 80, lcolor(gs8) lpattern(dash) msize(zero)) ///
			(pcarrowi 36 80 35 22, lcolor(gs8) lpattern(dash) msize(zero)) ///
			(scatteri 35 22 (9) "(22|35)" ,  mcolor(black) msize(medium) msymbol(O) mlabsize(small) mlabcolor(black)) ///
			(scatteri 45 17 (9) "(17|45)" , mcolor(black) msize(medium) msymbol(O) mlabsize(small) mlabcolor(black)) ///
			(scatteri 51 49 (2) "(49|51)" , mcolor(black) msize(medium) msymbol(O) mlabsize(small) mlabcolor(black)) ///
			(scatteri 45 82 (3) "(82|45)" , mcolor(black) msize(medium) msymbol(O) mlabsize(small) mlabcolor(black)) ///
			(scatteri 36 80 (3) "(80|36)" , mcolor(black) msize(medium) msymbol(O) mlabsize(small) mlabcolor(black)) ///
			(scatteri 35 22 (6) "A" ,  mcolor(black) msize(medium) msymbol(O) mlabsize(medlarge) mlabcolor(black)) ///
			(scatteri 45 17 (12) "B" , mcolor(black) msize(medium) msymbol(O) mlabsize(medlarge) mlabcolor(black)) ///
			(scatteri 51 49 (12) "C" , mcolor(black) msize(medium) msymbol(O) mlabsize(medlarge) mlabcolor(black)) ///
			(scatteri 45 82 (12) "D" , mcolor(black) msize(medium) msymbol(O) mlabsize(medlarge) mlabcolor(black)) ///
			(scatteri 36 80 (6) "E" , mcolor(black) msize(medium) msymbol(O) mlabsize(medlarge) mlabcolor(black)) ///
			(pcarrowi 20 30 34 34 (6) "Contour of Pareto Set", mlabcolor(black) lcolor(black) mcolor(black)) ///
			(pcarrowi 20 70 39 56 (6) "Uncovered Set", mlabcolor(black) lcolor(black) mcolor(black)), ///
				xtitle("") ///
				ytitle("") ///
				legend(off) ///
				yscale(range(0 100)) ylabel(0 (10) 100, labsize(vsmall) grid glwidth(thin)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, labsize(vsmall) grid glwidth(thin)) ///    
				graphregion(color(white) margin(medium)) ///
				plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
				note("Size of the Uncovered Set: 313 units" "Size of the Pareto Set: 780 units", size(vsmall)) ///
				xsize(4) ysize(4) ///
				name(Figure_2, replace)

*--------------------------------------------------------------------------------------------------------------------------
	
	*** Figure 3: Distributional properties of the policy space
	
		// Setting A
			twoway ///
				(contour t1_sd_points y x, scolor(white) ecolor(gs4) crule(linear) ccuts(0 (25) 500)) ///	
				(scatteri 51 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 47.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 50.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 46.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 51.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 45.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 52.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 42.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 55.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 41.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 56.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 40.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 58, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 60, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 60.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 60.2, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 59.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 58.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 58.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 57.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 41, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 56.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.75 55.9, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 43, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 45, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 46, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 47, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 48, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 50, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 51, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 52, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 53, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(pcarrowi 25 55 34 50 (6) "Contour of the uncovered set", color(black)) ///
				(scatteri 45 50, mcolor(black) msize(medium) msymbol(X)) ///
				(pcarrowi 70 70 46 51 (12) "Highest minimum payoff", color(black)) ///
				(scatteri 48 49, mcolor(black) msize(medium) msymbol(X)) ///
				(pcarrowi 80 30 49 48 (12) "Highest average payoff", color(black)) ///
				(scatteri 35 22 (6) "A", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 17 (12) "B", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 51 49 (12) "C", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 82 (12) "D", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 36 80 (6) "E", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)), ///
				zlabel(0 (50) 500, labsize(vsmall)) ///
				clegend( title("Standard" "Deviation" " of Payoffs" "in Tokens", size(vsmall)) width(4) forcesize region(margin(vsmall)) bmargin(tiny)) ///
				ztitle("") ///
				ytitle("") ///
				xtitle("") ///
				scheme(s1mono) ///
				yscale(range(0 100)) ylabel(0 (10) 100, nogrid labsize(vsmall)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, nogrid labsize(vsmall)) ///
				title("{bf: Setting A}", size(medsmall) color(black)) ///
				legend(off) ///
				aspect(1) ///
				graphregion(margin(vsmall)) ///
				plotregion(margin(vsmall)) ///
				name(Setting_A, replace)

		// Setting B
			twoway ///
				(contour t2_sd_points y x, scolor(white) ecolor(gs4) crule(linear) ccuts(0 (25) 500)) ///	
				(scatteri 51 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 47.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 50.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 46.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 51.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 45.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 52.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 42.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 55.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 41.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 56.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 40.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 58, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 60, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 60.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 60.2, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 59.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 58.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 58.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 57.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 41, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 56.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.75 55.9, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 43, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 45, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 46, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 47, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 48, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 50, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 51, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 52, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 53, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(pcarrowi 25 55 34 50 (6) "Contour of the uncovered set") ///
				(scatteri 42 68, mcolor(black) msize(medium) msymbol(X)) ///
				(pcarrowi 70 30 43 67 (12) "Highest minimum payoff", color(black)) ///
				(scatteri 40 71, mcolor(black) msize(medium) msymbol(X)) ///
				(pcarrowi 80 70 41 71 (12) "Highest average payoff", color(black)) ///
				(scatteri 35 22 (6) "A", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 17 (12) "B", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 51 49 (12) "C", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 82 (12) "D", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 36 80 (6) "E", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)), ///
				zlabel(0 (50) 500, labsize(vsmall)) ///
				clegend( title("Standard" "Deviation" " of Payoffs" "in Tokens", size(vsmall)) width(4) forcesize region(margin(vsmall)) bmargin(tiny)) ///
				ztitle("") ///
				ytitle("") ///
				xtitle("") ///
				scheme(s1mono) ///
				yscale(range(0 100)) ylabel(0 (10) 100, nogrid labsize(vsmall)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, nogrid labsize(vsmall)) ///
				title("{bf: Setting B}", size(medsmall) color(black)) ///
				legend(off) ///
				aspect(1) ///
				graphregion(margin(vsmall)) ///
				plotregion(margin(vsmall)) ///
				name(Setting_B, replace)

		// Setting C
			twoway ///
				(contour t3_sd_points y x, scolor(white) ecolor(gs4) crule(linear) ccuts(0 (25) 500)) ///	
				(scatteri 51 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 47.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 50.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 46.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 51.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 45.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 52.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 42.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 55.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 41.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 56.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 40.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 58, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 60, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 60.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 60.2, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 59.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 58.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 58.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 57.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 41, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 56.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.75 55.9, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 43, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 45, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 46, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 47, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 48, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 50, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 51, mcolor(blac#k) msize(vsmall) msymbol(O)) ///
				(scatteri 35 52, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 53, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(pcarrowi 25 55 34 50 (6) "Contour of the uncovered set") ///
				(scatteri 44 32, mcolor(black) msize(medium) msymbol(X)) ///
				(pcarrowi 70 70 45 33 (12) "Highest minimum payoff", color(black)) ///
				(scatteri 42 29, mcolor(black) msize(medium) msymbol(X)) ///
				(pcarrowi 80 30 43 29 (12) "Highest average payoff", color(black)) ///
				(scatteri 35 22 (6) "A", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 17 (12) "B", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 51 49 (12) "C", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 82 (12) "D", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 36 80 (6) "E", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)), ///
				zlabel(0 (50) 500, labsize(vsmall)) ///
				clegend( title("Standard" "Deviation" " of Payoffs" "in Tokens", size(vsmall)) width(4) forcesize region(margin(vsmall)) bmargin(tiny)) ///
				ztitle("") ///
				ytitle("") ///
				xtitle("") ///
				scheme(s1mono) ///
				yscale(range(0 100)) ylabel(0 (10) 100, nogrid labsize(vsmall)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, nogrid labsize(vsmall)) ///
				title("{bf: Setting C}", size(medsmall) color(black)) ///
				legend(off) ///
				aspect(1) ///
				graphregion(margin(vsmall)) ///
				plotregion(margin(vsmall)) ///
				name(Setting_C, replace)

			gen empty1 = _n/r(N)
			gen empty2 = _n/r(N)
			scatter empty1 empty2, name(empty, replace) msymbol(none)  /// 
				xlabel(,nogrid) ylabel(,nogrid) yscale(off) xscale(off) /// 
				graphregion(margin(zero)) ///
				plotregion(lwidth(none) margin(zero)) ///
				scheme(s1mono) ///
			  text(1 1 " " " " " " " " " " " " "  " " " " " " " " " " " " " "{it:Notes:}" "Setting A is used in Treatment 1, Period 1 & 2 |" "   Treatment 2, Period 1 | Treatment 3, Period 1." "Setting B is used in Treatment 2, Period 2." "Setting C is used in Treatment 3, Period 2." " " "The highest minumum payoff is the point from the policy space" "   that maximizes the payoff of the worst-off group member." "The highest average payoff is the point from the policy space" "   that maximizes the arithmetic mean of payoffs in a  group.",alignment(baseline) justification(left) size(small))
			  drop empty1 empty2
				
				
			graph combine  Setting_A Setting_B Setting_C empty, ///
				graphregion(color(white) margin(small)) ///
				cols(2) ///
				name(Figure_3, replace)
			graph display, ysize(10) xsize(11.5) scale(0.9) 

			graph close empty Setting_A Setting_B Setting_C 	
		serset clear

*--------------------------------------------------------------------------------------------------------------------------
		
		
	*** FIGURE B1: Representative payout function
	
		forvalues i=-100(10)100 {
		  local lab : display abs(`i')
		  local xlabs `"`xlabs' `i' "`lab'""'
			}
		display `"label string: [`xlabs']"'
		twoway ///
			(function y=(1000 - 5*abs(x)), range(-36 36) plotregion(margin(0)) color(black) clwidth(*1.5)) ///
			(function y=(9320*exp(-1* (abs(x) / 14.8))), range(-100 -36) plotregion(margin(0)) color(black) clwidth(*1.5)) ///
			(function y=(9320*exp(-1* (abs(x) / 14.8))), range(36 100) plotregion(margin(0)) color(black) clwidth(*1.5)) ///
			(function y=(820), range(0 36) plotregion(margin(0)) color(gs8) lpattern(dash) clwidth(*1.5)) ///
			(pcarrowi 950 25 950 15 (3) "Linear part", mlabcolor(black) lcolor(black) mcolor(black)) ///
			(pcarrowi 750 50 750 40 (3) "Exponential part", mlabcolor(black) lcolor(black) mcolor(black)), ///
			legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
			ytitle("{bf: Tokens}", size(medsmall)) ///
			xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(medsmall)) ///
			yscale(range(0 1000)) ylabel(0 (100) 1000, angle(0) labsize(small) grid glwidth(thin)) ///    
			xscale(range(-100 100) titlegap(*8)) xlabel(`xlabs', labsize(small) nogrid) ///   
			text(810 18 "{it: g}", place(s)) ///
			text(460 -40 "{bf:Payout function:}" " " ///
				"{it:Tokens} = 1000 - 5 *{it: distance}" "		if {it: distance} {&le} {it: g}" " " ///
				"{it:Tokens} = {it: k} * e{sup: -distance {&frasl} 14.8}" "		if {it: distance} {&gt} {it: g}" " " ///
				"{it:distance} = ((x - x{sup:opt}){sup:2} + (y - y{sup:opt}){sup:2}){sup:1/2}" " " ///
				"x{sup:opt}: x-coordinate of ideal point" ///
				"y{sup:opt}: y-coordinate of ideal point", ///
				place(se) box lcolor(black) fcolor(gs14) margin(small) just(left) size(small)) ///
			graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
			name(Figure_B1, replace)

*--------------------------------------------------------------------------------------------------------------------------


	*** FIGURE B2: Payout Functions

			forvalues i=-100(10)100 {
			  local lab : display abs(`i')
			  local xlabs `"`xlabs' `i' "`lab'""'
			}
			display `"label string: [`xlabs']"'

			* Setting A
			twoway ///
				(function y=(1000 - 5*abs(x)), range(-36 36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(-100 -36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(36 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///    
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player A}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_A_T1, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-36 36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(-100 -36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(36 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///  
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player B}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_B_T1, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-36 36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(-100 -36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(36 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///    
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player C}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_C_T1, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-36 36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(-100 -36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(36 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///  
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player D}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_D_T1, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-36 36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(-100 -36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(36 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///    
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player E}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_E_T1, replace)

			graph combine Player_A_T1 Player_B_T1 Player_C_T1 Player_D_T1 Player_E_T1, ///
				graphregion(color(white) margin(small)) ///
				cols(1) ///
				title(`"{bf: Setting A}"', size(medsmall) color(black)) ///
				note("Used in Treatment 1, Period 1 & 2" "Treatment 2, Period 1" "Treatment 3, Period 1", size(small)) ///
				name(T1_all, replace)

			* Setting B
			twoway ///
				(function y=(1000 - 5*abs(x)), range(-50 50) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(22000*exp(-1* (abs(x) / 14.8))), range(-100 -50) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(22000*exp(-1* (abs(x) / 14.8))), range(50 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///    
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player A}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_A_T2, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-55 55) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(29800*exp(-1* (abs(x) / 14.8))), range(-100 -55) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(29800*exp(-1* (abs(x) / 14.8))), range(55 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///   
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player B}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_B_T2, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-36 36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(-100 -36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(36 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///   
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player C}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_C_T2, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-12 12) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(2110*exp(-1* (abs(x) / 14.8))), range(-100 -12) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(2110*exp(-1* (abs(x) / 14.8))), range(12 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///  
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player D}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_D_T2, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-10 10) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(1865*exp(-1* (abs(x) / 14.8))), range(-100 -10) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(1865*exp(-1* (abs(x) / 14.8))), range(10 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player E}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_E_T2, replace)

			graph combine Player_A_T2 Player_B_T2 Player_C_T2 Player_D_T2 Player_E_T2, ///
				graphregion(color(white) margin(small)) ///
				cols(1) ///
				title(`"{bf: Setting B}"', size(medsmall) color(black)) ///
				note("     Used in Treatment 2, Period 2" " "  " ", size(small)) ///
				name(T2_all, replace)

			* Setting C
			twoway ///
				(function y=(1000 - 5*abs(x)), range(-10 10) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(1865*exp(-1* (abs(x) / 14.8))), range(-100 -10) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(1865*exp(-1* (abs(x) / 14.8))), range(10 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player A}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_A_T3, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-12 12) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(2110*exp(-1* (abs(x) / 14.8))), range(-100 -12) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(2110*exp(-1* (abs(x) / 14.8))), range(12 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player B}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_B_T3, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-36 36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(-100 -36) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(9320*exp(-1* (abs(x) / 14.8))), range(36 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player C}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_C_T3, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-55 55) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(29800*exp(-1* (abs(x) / 14.8))), range(-100 -55) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(29800*exp(-1* (abs(x) / 14.8))), range(55 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player D}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_D_T3, replace)

			twoway ///
				(function y=(1000 - 5*abs(x)), range(-50 50) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(22000*exp(-1* (abs(x) / 14.8))), range(-100 -50) plotregion(margin(0)) color(black) clwidth(*1.4)) ///
				(function y=(22000*exp(-1* (abs(x) / 14.8))), range(50 100) plotregion(margin(0)) color(black) clwidth(*1.4)), ///
				legend(off) xlabel(0 (0.1) 1, labsize(small)) ylabel(0(0.1)1, labsize(small) nogrid) ///
				ytitle("{bf: Tokens}", size(small)) ///
				xtitle("{bf: Euclidean Distance to Player's Ideal Point}", size(small)) ///
				yscale(range(0 1000)) ylabel(0 (200) 1000, angle(0) labsize(vsmall) grid glwidth(thin)) ///
				xscale(range(-100 100)) xlabel(`xlabs', labsize(vsmall) nogrid) ///   
				title("{bf: Player E}", color(black) size(medsmall)) ///
				graphregion(color(white) margin(medium)) xsize(4) ysize(3) ///
				name(Player_E_T3, replace)

			graph combine Player_A_T3 Player_B_T3 Player_C_T3 Player_D_T3 Player_E_T3, ///
				graphregion(color(white) margin(small)) ///
				cols(1) ///
				title(`"{bf: Setting C}"', size(medsmall) color(black)) ///
				note("     Used in Treatment 3, Period 2" " "  " ", size(small)) ///
				name(T3_all, replace)

			* ALL
			graph combine T1_all T2_all T3_all, ///
				graphregion(color(white) margin(small)) ///
				cols(3) ///
				name(Figure_B2, replace)
			graph display, ysize(12) xsize(12) scale(.55) 

			graph close Player_A_T1 Player_B_T1 Player_C_T1 Player_D_T1 Player_E_T1 T1_all Player_A_T2 Player_B_T2 Player_C_T2 Player_D_T2 Player_E_T2 T2_all Player_A_T3 Player_B_T3 Player_C_T3 Player_D_T3 Player_E_T3 T3_all
			serset clear
			
*--------------------------------------------------------------------------------------------------------------------------


	*** FIGURE B3: Average Payouts

		// Setting A 
			twoway ///
				(contour t1avg_points y x, scolor(white) ecolor(gs4) crule(linear) ccuts(0 (50) 1000)) ///	
				(scatteri 51 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 47.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 50.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 46.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 51.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 45.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 52.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 42.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 55.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 41.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 56.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 40.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 58, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 60, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 60.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 60.2, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 59.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 58.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 58.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 57.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 41, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 56.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.75 55.9, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 43, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 45, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 46, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 47, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 48, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 50, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 51, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 52, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 53, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(pcarrowi 25 55 34 50 (6) "Contour of the uncovered set") ///
				(scatteri 35 22 (6) "A", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 17 (12) "B", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 51 49 (12) "C", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 82 (12) "D", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 36 80 (6) "E", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)), ///
				zlabel(0 (100) 1000, labsize(vsmall)) ///
				clegend( title("Average" "Number" "of Tokens", size(vsmall)) width(4) forcesize region(margin(vsmall)) bmargin(tiny)) ///
				ztitle("") ///
				ytitle("") ///
				xtitle("") ///
				scheme(s1mono) ///
				yscale(range(0 100)) ylabel(0 (10) 100, nogrid labsize(vsmall)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, nogrid labsize(vsmall)) ///
				title("{bf: Setting A}", size(medsmall) color(black)) ///
				legend(off) ///
				aspect(1) ///
				graphregion(margin(vsmall)) ///
				plotregion(margin(vsmall)) ///
				note("Used in Treatment 1, Period 1 & 2 | Treatment 2, Period 1 | Treatment 3, Period 1", size(vsmall)) ///
				name(Setting_A, replace)

		// Setting B
			twoway ///
				(contour t2avg_points y x, scolor(white) ecolor(gs4) crule(linear) ccuts(0 (50) 1000)) ///		
				(scatteri 51 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 47.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 50.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 46.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 51.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 45.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 52.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 42.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 55.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 41.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 56.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 40.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 58, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 60, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 60.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 60.2, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 59.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 58.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 58.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 57.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 41, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 56.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.75 55.9, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 43, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 45, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 46, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 47, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 48, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 50, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 51, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 52, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 53, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(pcarrowi 25 55 34 50 (6) "Contour of the uncovered set") ///
				(scatteri 35 22 (6) "A", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 17 (12) "B", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 51 49 (12) "C", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 82 (12) "D", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 36 80 (6) "E", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)), ///
				zlabel(0 (100) 1000, labsize(vsmall)) ///
				clegend( title("Average" "Number" "of Tokens", size(vsmall)) width(4) forcesize region(margin(vsmall)) bmargin(tiny)) ///
				ztitle("") ///
				ytitle("") ///
				xtitle("") ///
				scheme(s1mono) ///
				yscale(range(0 100)) ylabel(0 (10) 100, nogrid labsize(vsmall)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, nogrid labsize(vsmall)) ///
				title("{bf: Setting B}", size(medsmall) color(black)) ///
				legend(off) ///
				aspect(1) ///
				graphregion(margin(vsmall)) ///
				plotregion(margin(vsmall)) ///
				note("Used in Treatment 2, Period 2", size(vsmall)) ///
				name(Setting_B, replace)

		// Setting C
			twoway ///
				(contour t3avg_points y x, scolor(white) ecolor(gs4) crule(linear) ccuts(0 (50) 1000)) ///	
				(scatteri 51 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 47.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 50.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 46.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 51.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 45.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 52.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 42.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 55.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 41.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 56.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 40.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 58, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 60, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 60.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 60.2, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 59.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 58.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 58.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 57.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 41, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 56.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.75 55.9, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 43, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 45, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 46, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 47, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 48, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 50, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 51, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 52, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 53, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(pcarrowi 25 55 34 50 (6) "Contour of the uncovered set") ///
				(scatteri 35 22 (6) "A", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 17 (12) "B", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 51 49 (12) "C", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 82 (12) "D", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 36 80 (6) "E", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)), ///
				zlabel(0 (100) 1000, labsize(vsmall)) ///
				clegend( title("Average" "Number" "of Tokens", size(vsmall)) width(4) forcesize region(margin(vsmall)) bmargin(tiny)) ///
				ztitle("") ///
				ytitle("") ///
				xtitle("") ///
				scheme(s1mono) ///
				yscale(range(0 100)) ylabel(0 (10) 100, nogrid labsize(vsmall)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, nogrid labsize(vsmall)) ///
				title("{bf: Setting C}", size(medsmall) color(black)) ///
				legend(off) ///
				aspect(1) ///
				graphregion(margin(vsmall)) ///
				plotregion(margin(vsmall)) ///
				note("Used in Treatment 3, Period 2", size(vsmall)) ///
				name(Setting_C, replace)


			graph combine  Setting_A Setting_B Setting_C, ///
				graphregion(color(white) margin(small)) ///
				cols(1) ///
				name(Figure_B3, replace)
			graph display, ysize(10) xsize(4.5) scale(1) 
			
			graph close Setting_A Setting_B Setting_C
			serset clear
				
*--------------------------------------------------------------------------------------------------------------------------


	*** FIGURE B4: Minimum payouts

		// Setting A 
			twoway ///
				(contour t1_min y x, scolor(white) ecolor(gs4) crule(linear) ccuts(0 (50) 1000)) ///	
				(scatteri 51 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 47.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 50.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 46.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 51.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 45.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 52.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 42.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 55.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 41.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 56.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 40.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 58, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 60, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 60.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 60.2, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 59.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 58.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 58.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 57.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 41, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 56.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.75 55.9, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 43, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 45, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 46, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 47, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 48, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 50, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 51, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 52, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 53, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(pcarrowi 25 55 34 50 (6) "Contour of the uncovered set") ///
				(scatteri 35 22 (6) "A", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 17 (12) "B", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 51 49 (12) "C", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 82 (12) "D", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 36 80 (6) "E", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)), ///
				zlabel(0 (100) 1000, labsize(vsmall)) ///
				clegend( title("Minimum" "Number" "of Tokens", size(vsmall)) width(4) forcesize region(margin(vsmall)) bmargin(tiny)) ///
				ztitle("") ///
				ytitle("") ///
				xtitle("") ///
				scheme(s1mono) ///
				yscale(range(0 100)) ylabel(0 (10) 100, nogrid labsize(vsmall)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, nogrid labsize(vsmall)) ///
				title("{bf: Setting A}", size(medsmall) color(black)) ///
				legend(off) ///
				aspect(1) ///
				graphregion(margin(vsmall)) ///
				plotregion(margin(vsmall)) ///
				note("Used in Treatment 1, Period 1 & 2 | Treatment 2, Period 1 | Treatment 3, Period 1", size(vsmall)) ///
				name(Setting_A, replace)

		// Setting B
			twoway ///
				(contour t2_min y x, scolor(white) ecolor(gs4) crule(linear) ccuts(0 (50) 1000)) ///		
				(scatteri 51 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 47.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 50.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 46.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 51.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 45.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 52.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 42.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 55.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 41.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 56.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 40.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 58, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 60, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 60.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 60.2, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 59.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 58.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 58.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 57.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 41, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 56.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.75 55.9, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 43, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 45, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 46, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 47, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 48, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 50, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 51, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 52, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 53, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(pcarrowi 25 55 34 50 (6) "Contour of the uncovered set") ///
				(scatteri 35 22 (6) "A", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 17 (12) "B", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 51 49 (12) "C", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 82 (12) "D", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 36 80 (6) "E", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)), ///
				zlabel(0 (100) 1000, labsize(vsmall)) ///
				clegend( title("Minimum" "Number" "of Tokens", size(vsmall)) width(4) forcesize region(margin(vsmall)) bmargin(tiny)) ///
				ztitle("") ///
				ytitle("") ///
				xtitle("") ///
				scheme(s1mono) ///
				yscale(range(0 100)) ylabel(0 (10) 100, nogrid labsize(vsmall)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, nogrid labsize(vsmall)) ///
				title("{bf: Setting B}", size(medsmall) color(black)) ///
				legend(off) ///
				aspect(1) ///
				graphregion(margin(vsmall)) ///
				plotregion(margin(vsmall)) ///
				note("Used in Treatment 2, Period 2", size(vsmall)) ///
				name(Setting_B, replace)

		// Setting C
			twoway ///
				(contour t3_min y x, scolor(white) ecolor(gs4) crule(linear) ccuts(0 (50) 1000)) ///	
				(scatteri 51 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 47.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.75 50.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 46.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.5 51.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 45.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50.25 52.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 50 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 42.75, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.75 55.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 41.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.5 56.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 40.25, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49.25 58, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 49 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48.5 60, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 48 60.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 47 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 46 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 45 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 44 60.2, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 37.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 43 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 42 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 41 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40.5 60.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 40 59.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 38, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39.5 59.5, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 39 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38.5 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 38 59.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 39, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37.5 58.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 37 58.3, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 40, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36.5 57.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 41, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 36 56.8, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.75 55.9, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35.5 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 42, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 43, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 44, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 45, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 46, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 47, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 48, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 49, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 50, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 51, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 52, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 53, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 54, mcolor(black) msize(vsmall) msymbol(O)) ///
				(scatteri 35 55, mcolor(black) msize(vsmall) msymbol(O)) ///
				(pcarrowi 25 55 34 50 (6) "Contour of the uncovered set") ///
				(scatteri 35 22 (6) "A", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 17 (12) "B", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 51 49 (12) "C", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 45 82 (12) "D", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)) ///
				(scatteri 36 80 (6) "E", mcolor(black) msize(medium) msymbol(o) mlabcolor(black)), ///
				zlabel(0 (100) 1000, labsize(vsmall)) ///
				clegend( title("Minimum" "Number" "of Tokens", size(vsmall)) width(4) forcesize region(margin(vsmall)) bmargin(tiny)) ///
				ztitle("") ///
				ytitle("") ///
				xtitle("") ///
				scheme(s1mono) ///
				yscale(range(0 100)) ylabel(0 (10) 100, nogrid labsize(vsmall)) ///    
				xscale(range(0 100)) xlabel(0 (10) 100, nogrid labsize(vsmall)) ///
				title("{bf: Setting C}", size(medsmall) color(black)) ///
				legend(off) ///
				aspect(1) ///
				graphregion(margin(vsmall)) ///
				plotregion(margin(vsmall)) ///
				note("Used in Treatment 3, Period 2", size(vsmall)) ///
				name(Setting_C, replace)

			graph combine  Setting_A Setting_B Setting_C, ///
				graphregion(color(white) margin(small)) ///
				cols(1) ///
				name(Figure_B4, replace)
			graph display, ysize(10) xsize(4.5) scale(1) 

			graph close Setting_A Setting_B Setting_C
			serset clear
